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Abstract 

We investigate the exclusive B — > ttI + £~ and B — * p£ + £~ decays in framework of 
the general two Higgs doublet model (model III), in which an extra phase angle in 
the charged-Higgs fermion coupling, i.e., a new source for CP violation exists. The 
CP violation for both decays are calculated and it is observed that the CP violating 
asymmetry in model III differs significantly than the one predicted by the standard 
model and model II which is a special case of model III. Furthermore, it is shown that 
the zero value of forward backward asymmetry Afb is shifted when compared with 
the SM value, which can also serve as the efficiency tool for establishing new physics. 



PACS numbers: 12.60. i, 13.20.-v, 13.25. Gv 



*e-mail: taliev@rorqual.cc.metu.edu.tr 
t e-mail : savci@rorqual . cc . metu.edu.tr 



1 Introduction 



Rare B meson decays, induced by flavor-changing neutral current (FCNC) b — > s(d) transi- 
tions, is one of the most promising research area in particle physics. Theoretical interest to 
the B meson decays lies in their role as a potential precision testing ground for the standard 
model (SM) at loop level. Experimentally, these decays will provide quantitative informa- 
tion about the Cabibbo-Kobayashi-Maskawa (CKM) elements V t d, V ts and Vtb- Besides 
these rare decay have the potential for establishing new physics beyond SM, such as two 
Higgs doublet model (2HDM), minimal supersymmetric extension of the SM (MSSM) |]J. 

Firstly the most reliable quantitative test of FCNC processes in B decays is expected to 
be measured in inclusive channels. In particular, the decays B — > X s ^£ + £~ are important 
probes of the effective Hamiltonian governing the FCNC transition b — > s(d)£ + £~ . The hope 
that B — > X s £ + £~ decay will be measurable in experiments in the near future, encourage 
extensive investigation of this process in the SM, 2HDM and MSSM P|-[15|. The matrix 



element of the b — > s£ + £~ contains terms describing the virtual effects induced by tt, cc and 
uu loops which are proportional to V t bV t * s , Vt, c V* s and Vb u V* u , respectively. Using unitarity 
of the CKM matrix and neglecting V& u V* u i n comparison to V t b V t * and V& c V* s , it is obvious 
that the matrix element for the b — > s£ + £~~ involves only one independent CKM factor 
V t b V t * s so that CP-violation in this channel is strongly suppressed in the SM. 

The situation is totally different for the b — > d£ + £~ decay, since all three CKM factors 
V t b Vtj, V c b V* d and V u b V* d , are all of the same order (in SM) and therefore can induce 
considerable CP-violating difference between the decay rates of the reactions b — > d£ + £~ 
and 6 -> d£+£~. 

It should be noted here that in presence of a much stronger decay b — > s£ + £~, the 
detection of the b — *> d£ + £~ decay seems to be more problematic. For this reason, in search 
of CP violation the corresponding exclusive decay channels B — > n£ + £r and B — > p£ + £~ 
are more preferable. In general, the inclusive decays are rather difficult to measure in 
comparison to the exclusive ones. CP-violating effects in inclusive b — > d£ + £~ and exclusive 
B — > n£ + £r , B — > p£ + £~ channels were studied within the framework of the SM in |T5|, IB 



The aim of the present work is to derive quantitative predictions for the CP violation 
in the exclusive B — > n£ + £~ and B — > p £ + £~ decays, in context of the general two Higgs 
doublet model, in which a new source for CP violation is present (see below). 2HDM model 
is one of the simplest extension of the SM, which contains two complex Higgs doublets, 
while the SM contains only one. In general, in 2HDM the flavor changing neutral currents 



(FCNC) that appear at tree level, are avoided by imposing an ad hoc discrete symmetry |L8 



One possible approach to avoid these unwanted FCNC at tree level is to couple all fermions 
to only one of the above-mentioned Higgs doublets (model I). The other possibility is the 
coupling of the up and down quarks to the first and second Higgs doublets, with the vacuum 
expectation values V2 and V\, respectively (model II). Model II is more attractive since its 
Higgs sector coincides with the ones in the supersymmetric model. In this model there 
exist five physical Higgs fields: neutral scalars H°, h°, neutral pseudoscalar A and charged 
Higgs bosons H ^ . The interaction vertex of fermions with Higgs fields depends on tan(3 = 
V2/V1, which is the free parameter of the model. The new experimental results of CLEO 
and ALEPH Collaborations [IS, 2C] on the branching ratio b — > sj decay impose strict 
restrictions on the charged Higgs boson mass and tan/3. Recently, the lower bound on these 
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parameters were determined from the analysis of the b — > 57 decay, including NLO QCD 
corrections [21], |2~2 |. Other indirect bound on the ratio m H ± /tan/3 come from i? — > Dtv t 
decay, where ra#± > 2.2tanf3 GeV |23|], and from the r lepton decays > 1.5 tanf3 GeV 
|2j. The consequence of an analysis without discrete symmetry has been investigated in 
a more general model in 2HDM, namely, model III |25|, |2f|. In this model FCNC appears 
naturally at tree level. However, the FCNC's involving the first two generations are highly 
suppressed, as is observed in the low energy experiments, and those involving the the third 
generation is not as severely suppressed as the first two generations, which are restricted 
by the existing experimental results. 

In this work we assume that all tree level FCNC couplings are negligible. It should be 
noted however that, even with this assumption, the couplings of fermions to Higgs bosons 
may have a complex phase e ld . In other words, in this model there exists a new source of 
CP violation that is absent in the SM, model I and model II. The effects of such an extra 

The constraints on the phase angle 9 



phase in the b — > S7 decay were discussed in 
in the product X u Xbb of Higgs-fermion coupling (see below) imposed by the neutron electric 
dipole moment, B° — B° mixing. p parameter and is discussed in [28]. 

The paper is organized as follows: In Section 2 we present the necessary theoretical 
framework. The branching ratios, CP-violating effects in the partial widths and forward- 
backward asymmetry for the above-mentioned exclusive decay channels are studied in sec- 
tion 3. Section 4 is devoted to the numerical analysis and concluding remarks. 



2 Theoretical framework 

Before presenting the necessary theoretical background, let us go through the main essential 
points of the general Higgs doublet model (model III). In this model, both Higgs doublets 
can couple to up and down quarks. Without loss of generality we can work in a basis such 
that the first doublet generates all the fermion and gauge boson masses, whose vacuum 
expectation values are 



( \ 



V y/2 ) 



2) ■ 



In this basis the first doublet <\>\ is the same as in the SM, and all new Higgs bosons result 
from the second doublet fa, which can be written in the following form 




V2 



V2H+ 
X ° 2 + iA° 



where G + and G° are the Goldstone bosons. The neutral Xi an d X2 are n °t the physical 
mass eigenstate, but their linear combinations give the neutral H° and h° Higgs bosons: 



Xi 
X2 



" - H°cosa — h°sina 



" - H°sina + h°cosa . 
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The general Yukawa Lagrangian can be written as 



(1) 



ta 2 (p, r\lf and , in general, are the non- 



where z, j are the generation indices, (f) = icr 2 0, r]^' 

diagonal coupling matrices, L = (1— 7 5 )/2 and R = (l+7 5 )/2 are the left- and right-handed 
projection operators. In Eq. (1) all states are weak states, that can be transformed to the 
mass eigenstates by rotation. After performing this rotation on the Yukawa Lagrangian, 
we get 



Y 



-H + U 



Vckm^R — i U VckmL 



v, 



(2) 



where U{V) represents the mass eigenstates of u, c, t (d, s, b) quarks. In the present 
analysis, we will use a simple ansatz for £, U+,T> 



V2m 



w 



(3) 



Also it assumed that is complex, i.e., Ay = |Ay| e ld , and for simplicity we choose ^ u ' v to 
be diagonal to suppress all tree level FCNC couplings, and as a result Ajj are also diagonal 
but remain complex. Note that the results for model I and model II can be obtained from 
model III by the following substitutions: 



X u = cot/3 Xbf, = — cot(3 for model I , 
X u = cot (3 Xbb = + tan/3 for model II , 



(4) 



and 9 = 0. 

After this brief introduction about the general Higgs doublet model, let us return our 
attention to the b —>■ d£ + £~ decay. The powerful framework into which the perturbative 
QCD corrections to the physical decay amplitude incorporated in a systematic way, is 
the effective Hamiltonian method. In this approach, the heavy degrees of freedom in the 
present case, i.e., t quark, W ± , H ± , h°, H° are integrated out. The procedure is to match 
the full theory with the effective theory at high scale /x = raw, and then calculate the 
Wilson coefficients at lower // ~ 0(m&) using the renormalization group equations. In our 
calculations we choose the higher scale as \i = mw, since the charged Higgs boson is heavy 
enough (m^± > 210 GeV see ||21|| ) to neglect the evolution from m#± to mw- 

In the version of the 2HDM we consider in this work, the charged Higgs boson exchange 
diagrams do not produce new operators and the operator basis is the same as the one used 
for the b — > d£ + £~ decay in the SM. For this reason in the model under consideration, the 
charged Higgs boson contributions to leading order change only the value of the Wilson 
coefficients at mw scale, i.e., 



<3 



2HDM i 



lm w 



(-O.B.DM ( \ 



rilHDM 

Wo 



Cf M (m^)+Cf [m w ) 
Cl M \m w ) + Cf ± (m w ) 



{m w ) = C b lQ M {m w )+C^ {m w ) 
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The coefficients C 2HDM {m w ) to the leading order are given by 

° 7 [mw) ~ X 24(x-l)3 + A{x-iy lUX 

y{7-hy-%y 2 ) , y 2 (3y - 2) 



72(y - 1)3 12(y - I] 
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, , / 2/(3 - 5y) y(3y - 2) \ 



+ 



Sm 2 #w svn-uw 

-19a; 3 + 25a; 2 -3a; 4 + 30a; 3 - 54a; 2 + 32a; - 8 4 
360r-l)3 + 18(x - l) 4 enx+ g 



B(m w ) + 



1 — Asin 2 9 w 
sin 2 9\Y 



C{m w ) 



+ I A« 



36(x-l)3 

1 — Asin 2 9 w xy 



sin 2 9 w 5 yy — 1 — jj* 
47y 2 - 79y + 38 3y 3 -6y 3 + 4 V 
108(y -l) 3 18(y-l) 4 



18(x- l) 4 
1 1 

TJTi: 



- v 



■ £ny 



(6) 



C M K) = —^-(B{m w )-C{m w )) 



where 



+ + , (7) 

stn 2 9 w 8 \ y-\ (y - 1) 2 / 



£(ar) = - 77 ^ rr + X inx , 



A(x-l) 4(x-l) 5 

c(x) = - x -( r\ + ^ +2 > 

V ; 4V3(a;-l) 2(x - l) 2 , 



m t 2 



2 ' 



and sin 2 9 w = 0.23 is the Weinberg angle. It follows from Eqs. (5-7) that among all 
the Wilson coefficients, only C 7 involves the new phase angle 9. We have neglected the 
neutral Higgs boson exchange diagram contributions, since Higgs boson-fermion interaction 
is proportional to the lepton mass. 

The effective Hamiltonian for the b — > d£ + £~ decay is [29-32] 



W = -4 Gf 
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|E<W<W + *«E<W [Oik) - on»)} 
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where 



A,, 



V, lh V. 



ubV ud 



v tb v; d 

and Ci are the Wilson coefficients. The explicit form of all operators Oi can be found in 
[29-32]. 

The evolution of the Wilson coefficients from the higher scale /i = mw down to the low 
energy scale fi = m;, is described by the renormalization group equation 



— C 



e//(A») 



q?" WOO , 



eff/ 



where 7 is the anomalous dimension matrix. The coefficient Cj (p) at the scale 0(m&) in 
NLO is calculated in [§1], f22|: 



C e 7 ff (m b ) = C° 7 (m b ) + ^^-C l 7 ^{m b ) 



47T 



where C^mt) is the leading order (LO) term and (nib) describes the NLO terms, 

whose explicit forms can be found in plf . In our case, the expressions for these coefficients 
can be obtained from the results of [21] by making the following replacements: 



\Y\ 



|A«| 



and XY* 



|AttAbJ e 



In the SM, the QCD corrected Wilson coefficient Cg(rrib), which enters to the decay am- 
plitude up to the next leading order has been calculated in [29-32] . The Wilson coefficient 
C10 does not receive any new corrections at all, i.e., C w {mi,) = C 2 Q DM (mw)- As we have 
already noted, in the version of the 2HDM we consider in this work, there does not appear 
any new operator other than those that exist in the SM, therefore it is enough to make 



the replacement C, 



SM, 



m w 



Cl HDM (m w ) in [29-32], in order to calculate C% HDM at m b 



scale. Hence, including the NLO QCD corrections, Cg(mb) can be written as: 

cm = cl HDM {») 



1 + ^u(s) 



7T 



+ g(m c , s) 3d(/x) + C 2 (/x) + 3C 3 (/x) + C 4 (/x) + 3C 5 (m) + C 6 (//) 



g(m c , s) - g(0, s)\ (3d(/i) + C 2 (ji)) - -g(0, s) (C 3 (/i) + 3C7 4 (/x)) 



- l -g (1, s) (4C 3 + 4C 4 + 3C 5 + C 6 ) - ^ (0, s) (C 3 + 3C 4 ) 
+ ^(3C 3 + C 4 + 3C 5 + C 6 ) , 



(9) 



where m c = m c /mb , s = p 2 /m|, and 



w(s) 



4 



-7T 



9 3 
5 + 4s 

3 (1 + 2s) 



Li 2 (s 

fed 



^n(s) - s) 



s - 



2s (1 + s) (1 - 2s) n ,^ 5 + 9s-6§ 2 
in (s) + 



3(l-sT(l + 2s) 



3(1 -s) (1 + 2s) 



(10) 
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represents the O (a s ) correction from the one gluon exchange in the matrix element of Og, 
while the function g (m c , s) arises from one loop contributions of the four-quark operators 
0\-Oq, whose form is 

8„ / 8 4 2 . 

-gM™0 + 27 + g^-g(2 + B) 



+ ^\l-y l \{Q{l-y l )[in + \^ = - i t^J + 9 ( y< - 1) 2arctan ^J_ | ,(11) 



where ?/j = 4mf/p 2 . 

The Wilson coefficients Cg receives also long distance contributions, which have their 
origin in the real uu, dd and cc intermediate states, i.e., p, u and J/ip, ip', ■ ■ ■■ In the case 
of the J / ip family this usually accomplished by introducing a Breit-Wigner distribution for 
the resonance through the replacement ([4-7,33]) 

( . ~ ( . ~ 3tt _ m Vi T(y t -> i+r) 
g (m c , s) ^ g {m c , s) — « 2^ T~2 2Y~- ^~ > 12 

where the phenomenological parameter k = 2.3 is chosen in order to reproduce correctly 
the experimental value of the branching ratio (see for example ||L6|| ) 

B[B -> J/^X -> ) = B(B -> J/^X) B( J/V» -> X£+r) . 

In order to avoid the double counting, in this work, as an alternative to the functions 
g (rh u , s) and g (m c , s) that describe the effects of uu and cc loops, we have used a dif- 
ferent procedure, in which these functions are expressed through the normalized vacuum 
polarization n^ ad (s) that is related to the experimentally measurable quantity 

D (T io t(e + e- -> hadrons) 

Rhad(s) = -— — — — - , 13 

cr(e+e — > p + p ) 

via the dispersion relation ( see |T5], [T7| for more detail). In this way it is possible to include 
the p, lu, J/ip, ip', ■ ■ ■ resonances into the differential cross section in an approximate way, 
consistent with the idea of global duality. In this approach the u and J / ip family resonances 
are well described through the Breit-Wigner form and p resonance is introduced by 

1 / * 2 \ 3 / 2 

iC= 4(1-4^) l^(^)! 2 . (14) 

where F w (s) is the pion form factor that is represented by a modified Gounaris-Sakurai 
formula (see [|4], [35 



The effective short-distance Hamiltonian for b — > d£ + £ decay [29-32] leads to the QCD 
corrected matrix element (when the d quark mass is neglected) 

m = |^^{c 9 e// J 7 ^ 

_ 2C 7 7 ^-dia^(l + l5 )b£^£ , (15) 

where p 2 is the invariant dilepton mass. In Eq. (12) all Wilson coefficients are evaluated 
at the p = nib scale. 
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3 The exclusive B — > n£ + £ and B —> p£ + £ decays 



In this section, we proceed to calculate the branching ratio and CP violating asymmetry 
in the B — > tt£ + £~ and B — > p£ + £~ decays. It follows from the matrix element of the 
b — > di + i~ that in order to be able to calculate the matrix element of the exclusive decay 
B -> M£ + £~, the matrix elements (M\d^(l + ^ 5 )b\B) and (M\dia, u p u (l + -f 5 )b\B) (M = n 
or p) have to be calculated. These matrix elements can be parametrized in the following 
way: 



(Ptt) <fy»(l - 75)6 B (p B ) 
n (p n ) di cr, v p v (l + 75)6 B (p B ) 

P(p p ,e) <fy*(l -75)6 



/ + (p 2 )(PB+^v + r(p 2 K, 

/ \ 2 / 2 o M JtIp 2 

(j>B + Pn),P -Vp. [ m B ~ m 7r) 



(16) 
(17) 



p(p,e) dia lu/ p v (l + 'y 5 )b B(p B 



- -e a c^V 2V{p2) 

— tpuXat P p Pb . 

F m B + m p 
~ ie *p( m B + m p )Ai(p 2 ) + i(pB+P P )p(e*p) 

+ wA e *p)^f \Mp 2 ) - Mp' 

= 4e^ ACT 5*Vp a T 1 (p 2 ) 



A2 (P 2 ) 
m B + rrip 

(18) 



+ 2% e*Jm 2 B - m 2 ) - (p B + p p ) fl (e*p) T 2 (p 2 



+ 2i(e*p) 



Pp - (Pb + P P )p 



m. 



Up 2 ) . (19) 



In all of the matrix elements above, p = Pb — Pm (M = n or p) and e* is the four- 
polarization vector of the p meson. Using Eqs. (15-19) we obtain for the matrix elements 
of the B -> irl+l~ and B -> p£ + £~ decays: 



M 



B — >7T 



Ga 
2^2% 



V tb V;A (2Ap^ + Bp,) l liX t + (2C Pn , + Dp,) h^i , 



(20) 



M B^p 



Ga 
2V2rt 



2A 1 e^ Xa e* u P y B + iB^l - iB 2 (e*p)(p B + P P )p- iB 3 (e*p)p, 



+ t^ld \2C x e, vXu e* v p x p p B + iD 1 e; - tD 2 (e*p)(p B + p P ), - tD 3 (s*p) Pp 



(21) 



where 



A = c e 9 ff r-c 7 



2rn b f T (p 2 ) 
m B + m-K 



B 

C 



m B — Tn n — p 
m B + m n 



= C 10 f 



+ 
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D = 

A l = 

B 1 = 

B 2 = 

B 3 = 

Ci = 

D x = 

D 2 = 



c 10 (f + + r) , 



m B + m p 



C e 9 ff {m B + m p )A l + AC 7 ^{m B - m 2 )T 2 , 



c «//_^ +4C ^ ( T2 + 



m B + mp 



m r — m A r 



-a" 2 -^(A 3 -A ) + AC 7 ^T 3 , 



y 9 o 
p2 



C 



V 



10 : , 

m B + m p 

Ci {m B + m p )Ai , 
A 2 



C 
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2m 



(22) 



^3 = C 10 ^(A 3 -A ) . 

Using Eqs. (20) and (21) and performing summation over final lepton and p meson polariza- 
tion (in the B — > pl + l~ case), we obtained the following results for the double differential 
decay rates (the masses of the leptons, in our case electron or muon, are neglected): 



G 2 o? \V tb V*\ 2 V\ 



dp 2 dz 2 11 7r 5 m B 



Xm 4 B (l - z 2 ) 



(\A\ 2 +\C\ 2 ) 



(23) 



dp 2 dz 



2 12 7r 5 m B 



2\m% 



+ 



x 



m%s(l + z 2 ) (|^| 2 + Id) 2 ) 
(A(l - z 2 ) + 8rs) (jBil 2 + |A| 2 ) - 2\m%{\ - r - s)(l - z 2 ) 
(Re (B^;) + Re (D^Z) ) + X 2 m%(l - z 2 ) ±- (\B 2 \ 2 + \D 2 



2r 



m% 



+8m B szVX[Re {B^C{) + Re (A x Dl) 



(24) 



where z = cosO , 9 is the angle between the three-momentum of the £ + lepton and that of 
the B meson in the center of mass frame of the lepton pair, A(1,tm, s) = 1 + r 2 M + s 2 — 



2r M - 2s - 2r M s, r M 



2 2 

and s = (M = 7i or p). The CP violating asymmetry 
m B m B 
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between B — > M£ + £ and B — > M£ + £ decays is denned as 

dT dT 

a _ dp 2 dp 2 

dp 2 <ip 2 



(25) 



where 



and 



dr (5 -> M£+r) 

dp 2 



rfr (b -> mi+i~) 

dp 2 dp 2 dp 2 

The differential decay widths B — > 7r£ + £~ and _B — > p can easily be obtained from Eqs. 
(23) and (24) by integrating over z. Finally we get the following results for CP violating 
asymmetry for the B — > n£ + £~ and B — > p£ + £~ decays 



A B c rip 2 ) 



(\Af + \Cf) 



+ f+h 



2m b 



m B + 



|/ + | 2 (ImA tt ) (Im^6) 



(Im^i) 7] 2 - (ImA u ) (Im£ 2 ) m + (ReA u ) (Im£ 2 ) t/ 2 



(26) 



W) ^ ^ { - 2 (ImA u ) (Im^ 2 ) 



where 



— Xm%s 



V 



TUB + Vtlr 



+ 



2\ 2 m% 
3r 



+ + 16rs^j (m B + m p ) |Ai| 2 - ^Am 



i (i z r-^ i4ii4> 



+ 



x 



+ 



2 (Im£i) r/ 2 - 2 (ImA u ) (Im£ 2 ) Vi + 2 (ReA u ) (Im&) t/ 2 
U\m%m b s TiV 8X 2 m%m b A 2 



3p 2 
2m 2 B m b / 4 



77iB + m p 



+ 



3rp 2 m B + m. 



T 2 + 



[m% - m 2 ) 



r f r p 



^-A + 16rs^j AiT 2 (m B + m p )(m 2 B - m 2 ) 



+ 



\\m B {\ -r - s)[ (m B + m p )^ aJt 2 + 
3 \ p l v \m B — m z p ) 

4m b (m B - m p ) 



T.) 



A 2 T 2 



p z 



(27) 



S p = ^ Xm %s(\A\ 2 + \C\ 2 )+^\ 2 m B (\B 2 \ 2 + \D 2 \ 2 ) 



1 

+ 27 



-A + 16rs) (IB^ + I^I 2 ) 
^Am|(l - r - s) ( (Re_Bi_B 2 ) + (Re£>i£> 2 ) ) 



(28) 
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In deriving these expressions, we have used the following parametrization 



Cg^ = £1 + A u £ 2 , 

Cf f = vi+im, (29) 

and assumed that all form factors are positive (see below). Interference of Cg^ and Cj^ 
terms gives new contribution to the CP violating asymmetry. The results for the CP 
asymmetry in model II can be obtained from Eqs. (26) and (27) by substituting Eq. (3) 
(i.e., r] 2 = 0). 

At the end of this section we present forward-backward asymmetry Afb, which involve 
different combination of the Wilson coefficients. The analysis of Afb is very useful in 
extracting precise information about the sign of the Wilson coefficients and the new physics. 
The forward-backward asymmetry is defined as 

i , dT f° r dT 
dz——— / dz- 



, / 2\ Jo dp 2 dz J-i dp 2 dz /onA 

= Tdz dT i f d P^~ ' 

Jo dp 2 dz J-i dp 2 dz 

The forward-backward asymmetry for the B — > n£ + £~ decay is zero, both in SM and 2HDM, 
in the limit m e — > 0. We can explain this fact briefly as follows. The hadronic current for 
B — > decay is a pure vector and the lepton current is also conserved when mi — > 0. 

The charge asymmetry (or Afb) is non-zero if there exist C-vilolating terms but such terms 
are clearly absent in the B — > n£ + £~. Using Eq. (24), the forward-backward asymmetry 
for the B — > p£ + £~ takes the following form: 



8m%sV\\ (ReBiC*) + (Re^D*) 



(31) 



FB - ' - - yj , 

Finally, we examine the CP-violating difference between A FB and A FB) i.e., 

5 Afb — Afb ~ Afb , 

with A F b being the forward-backward asymmetry in the anti particle channel, which can 
be obtained by the replacement 



Cg^(X u ) Cg^(X u — > A*) , 



whose explicit can easily be obtained from Eq. (31), with the above mentioned replacement 

-9 



of C e J f 



4 Numerical analysis 

Before presentation of our quantitative calculations and graphics, we would like to note that 
we have considered two different versions, namely model II and model III of the 2HDM, in 
our analysis. For the free parameters Ay, and X tt of model III, we have used the restrictions 
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coming from B — > X s 7 decay, B°-B° mixing, p parameter and neutron electric-dipole 
moment f28j, that yields |A&&| = 50, \X U \ < 0.03. 



The values of the main input parameteres, which appear in the expressions for the 
branching ratios, Afb and Aqp are: m& = 4.8 GeV, m c = 1.4 GeV, m T = 1.78 GeV, rriB = 
5.28 GeV,m n = 0.14 GeV. For B meson lifetime we take r{B) = 1.56" 12 s The 
values of the Wilson coefficients are, C x = -0.249, C 2 = 1.108, C 3 = 1.112 x HT 2 , C 4 = 
-2.569x 10~ 2 , C 5 = 7.4x 10~ 3 , C 6 = -3.144x 10~ 2 . Throughout the course of the numerical 
analysis, we have used the Wolfenstein parametrization of the CKM matrix elements, i.e., 



X u 



V ub V: d p(l-p)- V 2 + tr ] 



+ 0(X 2 



V tb V t * d (1 - p) 2 + n 2 

for which we have used the following three different sets of parameters, 

(0.3; 0.34) 

(P,V)=\ (-0.07; 0.34) 
(-0.3; 0.34) . 

Of course the explicit expressions for the form factors are needed in the present numerical 
analysis. In the current literature these form factors have been calculated in the framework 
of the three point QCD sum rule f3?J, relativistic quark model f38|, and light cone QCD 



sum rules [39-41]. In further numerical analysis we have used the light cone QCD sum 
rule predictions on the form factors. It should be noted that the light cone QCD sum rule 
predictions on the form factors are reliable in the region m 2 — p 2 ~ O (few GeV 2 ). In order 
to extend to the full physical region we have used best fitted expressions by extrapolating 
the numerical results with the condition that these approximate formulas reproduce the 
light cone QCD sum rule predictions to a good accuracy, in the above-mentioned region. 
The form of form factors which satisfy this condition can be written in terms of three 



parameters as 39, 40 



F(p 2 ) 



F(0) 



a F — + b F \ — T \ 



m. 



m 



B , 



where the values of parameters F(0), ap and bp for the relevant decays, B — ► 7r and B — * p, 
are listed in Table 1 (this Table is taken from [5S, 5C] 

Firstly we consider model II for numerical calculations. In Figs. (1) and (2) we present 
the dependence of the differential decay widths of the B — > 7re + e~ and B — > pe + e~ on 
p 2 for (p, rj) = (0.3; 0.34) at m H ± = 250 GeV and tan[3 = 1, with and without long 
distance contributions, correspondingly. In Figs. (3) and (4) we plot the variation of the 
CP- violating asymmetry Aqp with respect to p 2 , with the following set of parameters: 
(p, 77) = (0.3; 0.34) and tan/3 = 1. In both figures, the solid line corresponds to the SM 
case, dash-dotted and dotted lines represent the CP- violating asymmetry at two different 
values of the mass of charged Higgs boson m H ± = 250 GeV and 500 GeV, respectively. The 
total branching ratios for the B — > ire + e~ and (B — > pe + e~) decays at three different sets 
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of Wolfenstein parameters and at = 250 GeV are presented in Table 2. From Figs. 

1-4 we see that, in model II the dependences of the branching ratio and Acp asymmetry 
on p 2 are very similar to those predicted by the SM, but their magnitudes different in these 
models. These results are expected, since in model II, the charged Higgs contributions 
change only the values of the Wilson coefficients CV, Cg and C\g. In this version of the 
2HDM charged Higgs contributions give rise to constructive interference to the SM result. 
Therefore the branching ratio increases and CP asymmetry decreases. 

We presented in Table 2, the numerical values of the average values of the CP violating 

asymmetry (A C p), in the region 1 GeV 2 < p 2 < (mj/^ — 0.02 GeV^j , using the same 
values of the Wolfenstein parameters used in Figs. 3 and 4. 

The dependence of the forward-backward asymmetry on p 2 Afb{B — > pe + e~) is plotted 
in Fig. (5) for SM and model II, with the same set of parameters as in Fig. (1). It is 
observed that the value of p 2 at which A FB becomes zero is shifted in model II. Therefore, 
in future experiments, the determination of the value of p 2 at which A FB is zero can 
give unambiguous information about the presence of new physics. In Fig. (6) we plot 
the resulting difference in the forward-backward asymmetry for the values of Wolfenstein 
parameters (p, rj) = (—0.07; 0.2), with and without the long distance effects. From these 
figures we observe that, 5Fab for the B — > pe + e~ decay is positive in the non resonant 
region for all values of p 2 , both in SM and model II. 

Note that, the results we have presented for forward-backward asymmetry and its dif- 
ference are performed for (p, rj) = (—0.07; 0.2). However, for sake of completeness, we have 
gone through the same analysis for two different sets of the Wolfenstein parameters, namely, 
(p, rj) = (—0.3; 0.34) and (p, rj) = (—0.07; 0.34), as well as several different choices of tan/3. 
The numerical results and the relevant graphical presentations have demonstrated that, 
no remarkable differences have been observed among these different choices. In Figs. (7) 
and (8) we present the dependence of the CP asymmetry Aqp, integrated over p 2 , for the 
B — > 7re + e _ and B — > pe + e" decays on the phase angle 9 at mn± = 250 GeV, \Xbb\ = 50 
and |A tt | = 0.03, without the long distance effects in model III. From both figures, especially 
from B — > pe + e~ case, we observe that, the average CP asymmetry differs essentially from 
the one predicted by model II. In the region n/2 < 9 < 37r/2, the change in (A C p) is more 
than 2.5 times than that predicted by model II. This fact can be explained as, the charged 
Higgs and SM contributions interference destructively in the above-mentioned region of 9. 
It should be stressed that, depending on the value of the phase angle 9, the charged Higgs 
contributions can interfere with the SM results, either constructively or destructively. This 
case is absolutely different in model II, where the above-mentioned contributions interfere 
only constructively. The values of the branching ratios B — > tt£ + £~ and B — > p £ + £~ decays 
at different values of the phase angle 9 in model III are presented in Table 4. 

In conclusion, the exclusive B — > -n£ + £~ and B — > p£ + £~ decays are analyzed in the 
2HDM and it is found that, the CP violating asymmetry in model III differs essentially 
from the ones predicted by model II. 



12 



References 

[1] J. L. Hewett, In *Stanford 1993, Proceedings, Spin structure in high energy processes* 
463-475 (1993). 

[2] W. S. Hou, R. S. Willey and A. Soni, Phys. Rev. Lett. 58, 1608 (1987). 

[3] N. G. Deshpande and J. Trampetic, Phys. Rev. Lett. 60, 2583 (1988). 

[4] C. S. Lim, T. Morozumi and A. I. Sanda, Phys. Lett. B 218, 343 (1989). 

[5] B. Grinstein, M. J. Savage and M. B. Wise, Nucl. Phys. B 319, 271 (1989). 

[6] C. Dominguez, N. Paver and Riazuddin, Phys. Lett. B 214, 459 (1988). 

[7] N. G. Deshpande, J. Trampetic and K. Ponose, Phys. Rev. D 39, 1461 (1989). 

[8] P. J. O'Donnell and H. K. Tung, Phys. Rev. D 43, 2067 (1991). 

[9] N. Paver and Riazuddin, Phys. Rev. D 45, 978 (1992). 

[10] A. Ali, T. Mannel and T. Morozumi, Phys. Lett. B 273, 505 (1991). 

[11] A. Ali, G. F. Giudice and T. Mannel, Z. Phys. C 67, 417 (1995). 

[12] C. Greub, A. Ioannissian and D. Wyler, Phys. Lett. B 346, 145 (1995); 
D. Liu, Phys. Lett. B 346, 355 (1995); 
G. Burdman, Phys. Rev. D 52, 6400 (1995); 

Y. Okada, Y. Shimizu and M. Tanaka, Phys.Lett. B 405, 297 (1997). 
[13] A. J. Buras and M. Miinz, Phys. Rev. D 52, 186 (1995). 

[14] N. G. Deshpande, X. -G. He and J. Trampetic, Phys. Lett. B 367, 362 (1996). 

[15] S. Bertolini, F. Borzumati, A. Masiero, G. Ridolfi, 
Nucl.Phys. B 353, 591 (1991). 

[16] F. Kriiger and L. M. Sehgal Phys. Rev D 55,2799 (1997). 

[17] F. Kriiger and L. M. Sehgal Phys. Rev D 56, 5452 (1997). 

[18] S. Glashow and S. Weinberg, Phys. Rev D 15, 1958 (1977). 

[19] Talk by R. Briere, CLEO-CONF-98-17, ICHEP98-1011, in Proceedings of ICHEP98, 
Vancouver, Canada, July 1998; and in talk by J. Alexander, in Proceedings of 
ICHEP98, Vancouver, Canada, July 1998. 

[20] R. Barate et ai, ALEPH Collaboration, Phys. Lett. B 429, 169 (1998). 

[21] F. Borzumati and C. Greub, Phys. Rev D 58, 074004 (1998). 

[22] M. Ciuchini, G. Degrassi, P. Gambino, G.F. Giudice, Nucl.Phys B 527, 21 (1998). 

13 



[23] K. Kiers and A. Soni, Phys. Rev D 56, 5786 (1997). 
[24] A. Stahl, H. Voss, Z. Phys C 74, 73 (1997). 

[25] T.P Cheng and M. Sher, Phys. Rev D 35, 3484 (1987); ibid. D 44, 1461 (1991); 
W.S. Hou, Phys. Lett. B 296, 179 (1992); 

A. Antaramian, L. Hall, and A. Rasin, Phys. Rev. Lett 69, 1871 (1992); 
L. Hall and S. Weinberg, Phys. Rev D 48, 979 (1993); 
M.J. Savage, Phys. Lett B 266, 135 (1991). 



[26 
[27 
[28 
[29 
[30 
[31 
[32 

[33 

[34; 
[35 
[36 
[37 

[38 

[39 
[40 
[41 



D. Atwood, L. Reina, and A. Soni, Phys. Rev D 55, 3156 (1997). 

L. Wolfenstein and Y.L. Wu, Phys. Rev. Lett 73, 2809 (1994). 

D. Bowser-Chao, K. Cheung, W.-Y. Keung, |hep-ph/9811235| (1998). 

G. Buchalla, A. Buras, and M. Lautenbacher, Rev. Mod. Phys 68, 1125 (1996). 

A. J. Buras, M. Misiak, Miinz, and S. Pokorski, Nucl. Phys B 424, 374 (1994). 

B. Grinstein, R. Springer, and M. Wise, Nucl. Phys B 339, 269 (1990). 

M. Misiak, Nucl. Phys B393, 23 (1993), Erratum, ibid. B 439, 461 (1995); 
A. J. Buras and M. Miinz, Phys. Rev D 52, 186 (1995). 

A. I. Vainshtein, V. I. Zakharov, L. B. Okun and M. A. Shifman, 
Sov. J. Nucl. Phys 24, 427 (1976). 

F. Jegerlehner, Nucl.Phys.Proc.Suppl 51C, 131 (1996). 

T. Kinoshita, B. Nizic and Y. Okamata, Phys. Rev D 31, 2108 (1985). 

Particle Data Group, R. M. Barnett et ai, Phys. Rev D 54, 1 (1996). 

P. Colangelo, F. De Fazio, P. Santorelli and E. Scrimieri, 
Phys. Rev D 53, 3672 (1996). 

W. Jaus and D. Wyler, Phys. Rev D 41, 3405 (1990); 

D. Melikhov, N. Nikitin and S. Simula, Phys.Lett B 410, 290 (1997). 

P. Ball, J.High Energy Phys 9809:005 (1998). 

P. Ball, V. M. Braun, Phys.Rev D 58, 094016 (1998). 

T. M. Aliev, A. Ozpineci and M. Savci, Phys. Rev D 56, 4260 (1997). 



14 



Figure captions 



Fig. 1 Invariant mass squared (p 2 ) distribution of the branching ratio of the electron pair 
in the B — > 7re + e~ decay. Line 1 corresponds to the mass spectrum including the effects of 
p, uo and J/\& resonances, whereas line 2 corresponds to the non resonant invariant mass 
spectrum, in the SM. Analogously, lines 3 and 4 represents the same distributions, respec- 
tively, in the model II, at tanj3 — 1. In both models the Wolfenstein parameters are chosen 
to be (p,rp) = (0.3,0.34). 

Fig. 2 The same as in Fig. 1, but for the B — > pe + e~ decay. 

Fig. 3 CP-violating partial width asymmetry in the B — > 7re + e~ decay as a function 
of p 2 for the values of the Wolfenstein parameters (p, 77) = (0.3,0.34), including p, u; and 
J/^ resonances. Line 1 represents the SM. Lines 2 and 3 correspond to the model II case for 
the different choices of the charged Higgs boson mass m% = 500 GeV, and m% = 250 GeV, 
respectively. 

Fig. 4 The same as in Fig. 3, but for the B — > p e + e~ decay. 

Fig. 5 The dependence of the forward-backward asymmetry Aps on p 2 in the B — ■> pe + e~ 
decay. The Wolfenstein parameters are chosen to be (p, 77) = (—0.07,0.2). See Fig. 1 for 
the interpretation of the lines 1 to 4. 

Fig. 6 The CP-violating partial width asymmetry difference. 5fb = Afb — A-fb in 
the B — > peV decay for (p, 77) = (—0.07,0.2). See Fig. 1 for the interpretation of the 
lines 1 to 4. 

Fig. 7 The dependence of the CP violating asymmetry, integrated over p 2 , on the phase 
angle for the B — > 7re + e~ decay, in model III. In this figure the straight line corre- 
sponds to model II. The Wolfenstein parameters and the charged Higgs are chosen to be 
(p, 77) = (0.3,0.34) and m% = 250 GeV, respectively. 

Fig. 8 The same as in Fig. 7, but for the B — > pe + e~ decay. 
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F(0) 




b F 


A B C P 


0.26 ±0.04 


0.29 


-0.415 




0.22 ±0.03 


0.93 


-0.092 


yB^p 


0.34 ±0.05 


1.37 


0.315 


T B^p 


0.15 ±0.02 


1.41 


0.361 


rpB^p 


0.15 ±0.02 


0.28 


-0.500 


rpB^p 


0.10 ±0.02 


1.06 


-0.076 


f?-** 


0.30 ±0.04 


1.35 


0.270 


J-'B — >7T 

It 


-0.30 ±0.04 


1.34 


0.260 



Table 1: 





B(B -> vre+e-) 


B(B -> pe+e-) 


(p; v) 


SM 


THDM 


SM 


THDM 


(±0.3; 0.34) 


3.27 x 10~ 8 


4.11 x 10~ 8 


5.99 x 10- 8 


8.45 x 10~ 8 


(-0.3; 0.34) 


3.31 x 10" 8 


4.15 x 10" 8 


6.00 x 10~ 8 


8.46 x 10" 8 


(-0.07; 0.34) 


3.30 x 10~ 8 


4.14 x 10~ 8 


6.00 x 10~ 8 


8.46 x lO" 8 



Table 2: 
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m H ± = 250 GeV 


-0.048 


-0.030 


m H± = 500 GeV 


-0.054 


-0.031 



Table 3: 



6 


B(B -> vre+e-) 


B(B -> peV) 





3.12 x 10~ 8 


7.41 x 10~ 8 


tt/4 


3.16 x 10~ 8 


7.08 x 10~ 8 


vr/2 


3.26 x 10~ 8 


6.34 x 10~ 8 



Table 4: 
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